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Abstract 

A criterion is derived for the determination of parameter domains of first 
order phase transitions in quantum mechanical tunnehng models. The criterion 
is tested by application to various models, in particular to some which have been 
used recently to explore spin tunneling in macroscopic particles. In each case 
agreement is found with previously heuristically determined domains. 
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1 Introduction 

It is well known that at temperatures close to zero, decay rates of nietastable 
states are determined by quantum tunneling processes whose dynamics is de- 
scribed by classical configurations such as vacuum or periodic bounces, but that 
with increasing temperature thermal activation becomes more and more impor- 
tant, and beyond some critical or crossover temperature Tc becomes the decisive 
mechanism. It has been realised recently that this change or crossover can be 
regarded as a transition from quantum to classical behaviour which in turn can 
be looked at like a phase transition in e.g. the thermodynamics of gases. This 
is a fascinating topic which has attracted considerable interest recently and ap- 
plication particularly in the area of macroscopic spin systems. The general idea 
applies also to other cases where instead of a decay rate a shift in energy is 
involved, such as a level splitting in quantum mechanics which is described by 
vacuum or periodic instantons. 

The characteristic way in which phase transitions appear in quantum mechan- 
ical tunneling processes has been worked out in ref. |jl|. In particular, a sharp first 
order transition is there shown to appear in the plot of action S versus temper- 
ature T which is very analogous to that of the free enthalpy versus pressure of 
a van der Waals gas, whose equation of state plotted as pressure versus volume 
corresponds to the plot of the period P{E) (of the periodic bounce or instanton) 
versus energy in the consideration of quantum mechanical systems. It is well 
known that in the case of a periodic problem in Euclidean time the derivative 
of the action with respect to the energy E is the negative of the oscillation time 
t{E) or period P{E) at that energy which has to be identified with —h/T. If 
t{E) is a monotonically decreasing function with increasing E, one has a second 
order transition; if t{E) has a minimum and then rises again within the domain 

< E < barrier height, one has a first order transition, i.e. in this case there is an 
energy E^ within this range with T^ = h/T{Ec) at which the first order transition 
takes place. A criterion for a first order transition can be obtained by studying 
the Euclidean time period in the neighbourhood of the sphaleron configuration at 
the peak of the potential barrier, i.e. at the bottom of the well of the inverse po- 
tential, as advocated in ref. [^ . If the frequency of oscillation about the sphaleron 
point is ujg and oscillations different from Us about it are possible, a first order 
transition requires u"^ > uj"^ ot r — t^ < 0. Thus this is a condition one can use 
as a criterion for a first order transition. Explicit and solvable examples of both 
types of transitions in quantum mechanical models have been given recently in 
refs.p, ^, |, 0. It was shown there that simple models, such as those of the 

1 + 1 dimensional scalar field theories with the familiar double well or periodic 
or sine-Gordon potentials, do not lead to first order transitions. However, as 
shown in ref. [Q, this changes if the mass is allowed to become field dependent, 
a case which occurs, for instance, in spin tunneling theories, or if an externally 
applied field is introduced as in refs.0, ^, H, 0], as again in the context of spin 



tunneling phenomena. The introduction of such extra dependences imphes, of 
course, the introduction of appropriate parameters. Thus in ref. it was shown, 
for example, that in the case of a periodic potential and an effective mass with a 
field dependence multiplied by a parameter A, the period P{E) is monotonically 
decreasing for A < |, but for A > | it has a minimum and then rises again, and 
so leads to a first order transition. It is therefore of interest to find criteria for 
the determination of domains of such parameter values in which they lead to one 
type of phase transition or the other. A significant step forward in this direction 
on the basis of higher dimensional arguments, was achieved in ref.|@], where an 
expansion about the time independent saddle point configuration is considered. 
However, many considerations of basic interest are best illustrated by quantum 
mechanical models, which serve as useful prototypes for more complicated theo- 
ries, and it is natural to enquire about criteria for the occurrence of first order 
transitions in these, and hence to verify the conditions obtained heuristically as 
in the models referred to above 0, ^, ^. Thus in the following we consider a 
scalar theory with a field dependent mass and suitable nonlinear potentials. Ex- 
panding the static field (quantum mechanically a coordinate) about the constant 
sphaleron configuration (at the top of the barrier) and proceeding to second order 
perturbation theory, we obtain the deviation of the frequency from its sphaleron 
value (lowest order perturbation theory) and hence the criterion for a first or- 
der transition. We then consider various models considered earlier |^ ^, |^, |^, |^ 
and demonstrate the agreement of the prediction of the criterion with the results 
obtained there heuristically. We add that theories with effective field dependent 
masses are not unusual; typical examples are, for instance, Skyrme theories p. 
The instability of the sphaleron configuration at the top of the barrier in various 
familiar models has been shown, for instance, in 0. 

2 The condition for a first order transition 

We consider a Euclidean action integral of the form 
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S= I dr 
with Euler-Lagrange equation 



-M{q)q' + V{q) 



(1) 



,^, ,.. ldM{q) .2 dV{q) ,^, 

and g(r + P) = q{T), where P = j; is the period and T the temperature. Since 
the solutions near the sphaleron solution at the top of the barrier shown in Fig.l 
have information on the order of the "phase" transition between quantum and 
thermal activity regimes, we confine ourselves to this region in solution space. 



We define (as for a maximum of V{q) at qs) 

V"{qs 



-u, = 



M{qs 



(3) 



and set in eq.(0) 

(l = qs + riir) (4) 

wliere qg is the splialeron position with V'{qs) = 0. Expanding M{qs + ^(t)) and 
V{qs + i]{t)) and their derivatives in powers of i]{t), we obtain the fluctuation 
equation 



M(g,) + M'iq.Ur) + ^M'^q^Wir) 



riir) 



1 
+ 2 



1 



M'iqs) + M"iqMr) + -M'"(g,)r^2(r) + 



V\r) 



= V'iqsHr) + lv"\q,WiT) + lv""{qsW{r) + ■■■ 

We are interested in small fluctuations i]{t). The first order equation 

MiqsMr) - V"{qs)vir) = 

has the even solution 

?7(r) = acos(cjor), cUq = lu"^ 



(5) 



(6) 



(7) 



This is the solution we choose for the following reason, which we explain for 
simplicity for the case M{q) = 1. The classical equation 

.. _ dV 
dq 

is invariant under time translations r -^ r+ const. The operator of the fluctu- 
ation equation or second variational derivative of the action S at the classical 
configuration q^ must therefore possess a zero eigenvalue with eigenfunction -^. 
This operator is also obtained by differentiating the classical equation and implies 
the fluctuation equation 



l^-^"("))* 



^ii^i 



If qc is the constant g^, and if this is a saddle point, we must have (at and 
around g^) uf = 0,'?/'i odd (the first excited state), and Wg < 0, and ipo even (the 
ground state). Thus in this case ipo is proportional to cos(a;oT). This explains 
our ansatz (|^) which is the component of a general fluctuation in the direction 
of the negative eigenmode of the fluctuation equation. The fluctuation equation 
for the case M{q) ^^const. has been discussed in |T^ . 



For the solution in the next order of perturbation theory we set 

ri{T) = acos{ujT) + a^?7i(r), u'^ = uJq + aAiu'^ 



(8) 



We insert this again into eq.(^) and retain only terms up to and including those 
of Ola"^). Also we reexpress powers of cos(ct;r) in terms of cos{nuJT) where n is 
an integer and use eq.(|3). Then 



M{qs)^-V"{q, 



Vii^) 



+ a 



2 



-'M'M^v'M) ^ (s^msA + Z!M\ „,p,^,) 



(9) 



We now expand the fluctuation ?7i(r) in terms of lowest order functions, i.e. we 

set 

Vii^) = Yl CnCOs{nUJoT) 
n=0,l,2,-- 

and use the orthonormality of the latter, i.e. 



cos{muoT) cos{nu!oT)dT = —6mn, P = — 

2 ujo 



The remaining steps are standard in perturbation theory. Thus, multiplying 
eq-® by cos(u;ot) and integrating one obtains immediately that 



A,u' 







and the fluctuation rji^r) is obtained to be 

'7i(^) = 5'i + fi'2 cos(2u;oi") 
with 

u'M'{qs) + V'"{qs) , , 3M'{q,)u' + V'"{q, 



9i[^) 



4\/"(g. 



92[iO) 



4 [AM{q,)uji + V"{q,)] 



(10) 



(11) 



uJq. Thus, in order to obtain a 



where in these first order expressions tu^ 

nonvanishing deviation from the sphaleron value, we have to proceed to the next 

order of perturbation theory. Hence we set 



?7(r) = acos(co'r) + a^rji^r) + a?rj2{T), cj^ = cUq + aAitu^ + a^A^uj'^ 



(12) 



and insert this into eq.(^) and expand up to and including terms of 0{a?'). The 
procedure is the same as before but the calculations are now much more involved. 



so that we can only cite a main intermediate result. The equation corresponding 
to eq.(|]) above and multiplied by cos(a;or) and integrated over yields the equation 



-aM{qs)a^A2Uj' 



P 



^'"fe)(| + f)i^ + V'fe)^i^ 



-2M'iqs)co'g2j - ^M"iq,)co'^ 



(13) 



From this the condition of a first order phase transition, i.e. A2Ci;^ > 0, becomes 



V"\qs){gi + 



92 



-y'"'iQs 



+M'{qs)u'g2 + M'{qs)uj'ig, + ^ 



92 



+lM"iq,)u;' 



<0 



(14) 



UIQ 



The first two terms of the expression on the left can be seen to be contained in 
the appropriate formula of ref. (there eq.(31)) where M = const. 

3 Application of the criterion to various models 

In the following we focus our interest on the criterion for a first order transition 
in various model theories and not on the physical contexts of these theories. We 
therefore do not elaborate on the latter which are discussed in the appropriate 
sources cited. 

1. A spin tunneling model with variable mass but without applied magnetic 
field 

In ref.H a large-spin tunneling model was considered in which 



M(0) 



1 



V [^] = K2s{s + 1) sin^ (^ (15) 



2Ki (l- Asin^, 
Thus, identifying with our previous coordinate q, we have 

0s = |, V[^s]=K2s{s + l), V'[^s]=0, 



V" [0.1 



-2K2s{s + l), V" 



0, V" 



8K2s{s + 1) 



and 




















M [0,] = 


1 




M' [0,] = 0, M" 


4>s 


= - 


A 






2K,{1- 


A)' 


Kiil- 


A)^ 


From 


eq.d) we 


then obt 

a 


ain 

0^ = 


AKiK2s{s + l){l - 


-A) 








It follows that 


91 = 0,^2 


= 0, 


and eq.(14) becomes 










K2S{S 


+ 1)+ '- 

v 


A 


\ 1 

K K "i 


s + 


1)(1 


-A)< 


0, 




) ii - xy"'"''^ 


i.e. 








->l 











(16) 



(17) 

This is precisely the condition found heuristically in ref . [^ for the existence 
of a first order transition in the model. The result implies, of course, the 
necessity of a significant field dependence of the effective mass to generate 
the first order transition. In the next model this effect is produced by an 
applied magnetic field. 

2. A spin tunneling model with constant mass and applied magnetic field 

In ref. a large-spin tunneling case was considered in which (in our present 
notation) 

M = const., V [(p] = {s + -fD (hi sinh^ <p - 2h^ cosh 0) (18) 

Here h^ = i2s+i)d ' where H^ is the magnetic field applied in the direction 
of coordinte x which in our notation is also (p = x. In this case one finds 

0, = 0, V[<j)s] = -2{s + ^yDh,„ V^'[0.]=O, 

V" [0,] = -2(s + ^Yh^ii - K), V" [0,] = 0, 



and 

Then from eq. 



V""[ct>s] = 2{s + ]^fDh,{Ah^-l] 



M = M,M' = M" = 



^o = ^is + ^fh,{l-h,) (19) 



and gi = 0, g2 = and the inequality (JT^) becomes 

1 _ 1 



and hence 



^■2{s + ^)^Dh,{4h^~l)<0 



K. < -7 



(20) 



as found in ref.0, there expressed as H^ < sD/2. 



3. A spin tunneling model with variable mass and applied transverse magnetic 
field 

In ref.|^ another and more comphcated spin tunnehng model was investi- 
gated in which the magnetic field is applied in a direction transverse to that 
of the easy axis and the potential has the shape of a periodically recurring 
asymmetric twin barrier (i.e. a large one followed by a small one) with 



M{<t>) 

V[<j>] 



2K [1- X sin^ (h + aX sin 



KXs^ (sin0 — a) 



(21) 



In this case we have to distinguish between transitions in the large barrier 
and those in the small barrier. We therefore consider these separately. 

(i) Large harrier transitions 

In this case the sphaleron angle is (superscript L referring to the large 
barrier) 



and there 



V" 
and 



V 



-2KAs2(l + «), V" 
1 



-TT 

2 



KXs'^{l + af, V 

= 0, V"" 



2KXs'^{4: + a) 



MiVs 



2K[l-Xil + a)y ^'^^^) 



M' 



ll( aL\ 



X 



{2 + a) 



2K [1 - A (1 + a)X 



In this case we obtain with eq. 



UJ, 



^' - 4i^2^s'(l + «)[!- A (l + «) 



(22) 



Again gi = 0, g2 = and the condition for a first order phase transition is 
found to be 

^ ^ (a + 1) (3a + 8) ^^^^ 

The parameter domains for first order and second order transitions through 
the large barrier are shown in Fig. 2. 

(ii) Small barrier transitions 

In the case of transitions through the small barrier we have 0; 



I and 



V 



and 



KXs^ (1 



V" 



M{(t>l 



a] 



v 



0,1/" 



2K[l-X{l-a 



0, V\^ 

= 2K\s^ (4 - 



-2K\s^ (I 



a] 



a] 



M' 



i'(aS\ 



X 



(2-«) 



2K[l-X{l-a)Y 



From eq. (^ we find 



UJ, 



S2 



AK^Xs'^ (1 - a) [1 - A (1 - a)] 



(24) 



Again Qi = 0,g2 = and the condition for a first order transition through 
the small barrier is found to be 



A> 



a 



8 - 3a 1 



a] 



(25) 



The parameter domains for first and second order transitions through the 
small barrier are shown in Fig. 3. The parameter domains for transitions 
through either barrier are shown in Fig. 4. 

4. A spin tunneling model with variable mass and applied longitudinal field 

In ref . [0 yet another model was investigated in which the applied magnetic 
field assumes a longitudinal, i.e. parallel, direction with respect to the easy 
axis. In this case it is also possible to find explicit classical configurations 



m, but we do not consider these here. In this model 

M(0) ^ 

V[0] 



2Ki (l - A sin^ - ^ cos 
K2s'^ f sin^ (J) + a cos + a 



(26) 



In this case the sphaleron value of is 0s = arccos f and 



V 



a 



K2S^ ! + « + —, V 



V" 



-Kos' 2 - 



a 



V" 



-3K2S^a\ 1 



a" 



V" 



K2s' ( 8 - -«^ 



and 



M(0,) 



M" 



M'i 



A aA/l - 



A (K-2)(l-A) + ^(l-^) 



2Ki 



[i-xy 



In this case 



ul = 2KiK2S^{l-\) 2 



a 



The expressions for gi and g2 are now found to be 



a 



9i 



92 






a 



1 02 



A 



1-A 



1 



a 
T 



(27) 



(28) 



The condition for a first order phase transition is now more complicated in 
that it reduces to the inequality 

>,2\ 



where 



Setting 



re'''^'-^H'^Y}<° 



1-A I 4 



(29) 



(i 



3a' 



1± 



\ 



3 o / "' 
1 - -a2 1 + — 

4 I 2 , 



the inequality 



is then ^- < ^ < ^+, which reduces to 



e-+ 1-f 



< A< 



^4 



1-T 



(30) 



In Fig. 5 we plot the boundary separating the parameter domain of first 
order transitions from that of second order transitions. 
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4 Conclusions 

In the above we have derived a criterion for the parameter domains of first order 
phase transitions for a wide class of quantum mechanical tunneling models. By 
application to specific models which received attention recently, particularly in 
the discussion of the tunneling of large or macroscopic spins, we demonstrated the 
agreement of the prediction of the criterion with previous heuristically obtained 
results. We envisage that considerations like those presented above and in ref.p| 
can be useful in wider contexts, and also in field theory models, such as those 
with Skyrme terms which imply an effective field dependent mass. 
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Figure captions 



Fig.l 

The potential barrier with sphaleron coordinate Qs 



Fig. 2 

The parameter domains for first and second order transitions through the 
large barrier in the spin tunneling model with field dependent mass and field 
transverse to the easy axis 



Fig. 3 

The parameter domains for first and second order transitions through the 
small barrier in the spin tunneling model with field dependent mass and field 
transverse to the easy axis 



Fig. 4 

The parameter domains for first and second order transitions in the model 
with field dependent mass and transverse magnetic field 



Fig. 5 

The parameter domains for first and second order phase transitions in the 
model with applied longitudinal field 
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Fig. 5 
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